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We consider the analytical properties of the eigenspectrum generated by a class of central poten- 
tials given by V{r) = —a/r + br^, 6 > 0. In particular, scaling, monotonicity, and energy bounds are 
discussed. The potential V{r) is considered both in all space, and under the condition of spherical 
confinement inside an impenetrable spherical boundary of radius R. With the aid of the asymptotic 
iteration method, several exact analytic results are obtained which exhibit the parametric depen- 
dence of energy on a, b, and R, under certain constraints. More general spectral characteristics are 
identified by use of a combination of analytical properties and accurate numerical calculations of the 
energies, obtained by both the generalized pseudo-spectral method, and the asymptotic iteration 
method. The experimental significance of the results for both the free and confined potential V{r) 
cases are discussed. 
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I. INTRODUCTION 



The model for a hydrogen atom, HA, confined in an impenetrable sphere of finite radius R was originally introduced 
[1] to simulate the effect of high pressure on atomic static dipole polarizability. Sommerfeld and Welker [2] formulated 
the wave function solutions for this potential in terms of confluent hypergeometric functions, and underlined the 
application of this model for the prediction of the line spectrum originating from atomic hydrogen in the outer 
atmosphere. An algorithm for obtaining nearly exact energy calculations for a spherically confined hydrogen atom 
has been published 3]. On the other hand, regular soft confinement of the Coulombic systems has been developed 
by superimposing Debye screening [?j . Such a confining potential has been successful in explaining [21 [B] the shift in 
frequency of the x-ray spectral lines emitted by laser-imploded plasmas in the limit of high plasma density, whereby 
the effective potential assumes the form given by the Coulomb plus oscillator potential. The harmonic potential can 
be considered here as giving rise to the confinement of the Coulomb system with soft boundary walls. A variety of 
other model potentials leading to the confinement of electrons in atoms and molecules have been proposed, in order 
to explain the behavior of the novel artificial nanostructures, such as quantum wires and quantum dots, atoms and 
molecules embedded inside fullerenes, zeolites and liquid helium droplets, and, in addition, to simulate the interior of 
a giant planet. A comprehensive review covering of the development and applications of confining model potentials 
has been recently published [71 |S] . Under the confinement effect of an impenetrable spherical cavity of radius R, the 
hydrogen atom and the isotropic harmonic oscillator, IHO, potentials have been studied, independently, and their 
spectral characteristics have been analyzed [21 [10] in terms of useful quasi exact results. In the following text, we shall 
denote the spherically confined hydrogen atom as SCHA and the spherically confined isotropic harmonic oscillator 
as SCIHO: in both cases, the eigenstates are labelled as (i/, ^), — 1, 2, 3, • • • ,£ — 0, 1, 2 • • • , in terms of which the 
number of radial nodes for a given £ becomes v — £ — 1. 

For the free Coulomb plus oscillator potential, a few exploratory calculations have been reported earlier [Tmi4j . In 
view of the importance of the conjoined Coulomb and harmonic oscillator potential, it useful to study the general 
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behavior of this potential under the confinement due to an impenetrable spherical cavity, as a function of the radius 
i?, where the free state is represented by R — >oo. In this paper we consider a general spherically symmetric model 
of atomic system confined by (i) the presence of a harmonic-oscillator potential term and in a representative set of 
cases also (ii) containment inside an impenetrable spherical box of radius R. In atomic units h = m = e = 1 the 
Hamiltonian for the model system is given by 

H = -lA + V{r), V{r) = -'^+br\ (1) 
2 r 

where a and b are coupling parameters. We shall always assume that b > and for the most part, we shall assume 
that the Coulomb term is also attractive, a > 0; we shall also consider the repulsive case a < 0, in which H becomes 
a model, for example, for a system composed of a pair of confined electrons. 

We shall now present a brief review of the known results defining the spectral characteristics of the two confined 
systems SCHA and SCIHO. It is well known that the so called accidental degeneracy of free HA is removed in the 
SCHA. As i? — > 0, the energy levels £{1/, £) increase in magnitude such that the higher £ states get relatively less 
destabilized. There exists a critical value of R above which E{v^t) > 0. Further, two additional kinds of degeneracies 
arise [5]. They result from the specific choice of the radius of confinement R, chosen exactly at the radial nodes 
corresponding to the free HA wave functions. In the incidental degeneracy case, a given confined {v = £ + 1,1) state 
becomes iso-energic with {1/ — £ + 2,£) state of the free HA with energy —l/{2{£ + 2)^} a.u., at the same R. In the 
simultaneous degeneracy case, on the other hand, a certain pair of confined states at the common radius of confinement 
R that is prescribed in terms of the location of the radial node in a specific free state of HA, become iso-energic. For 
example, for all > ^ + 2, each (i^, ^) SCHA state is degenerate with {i^ + 1,£ + 2) state, when both of them are 
confined at R — {£ + !){£ + 2), which defines the radial node in the free {£ + 2, £) state. Both these degeneracies have 
been shown 9 to result from the Gauss relationship applied at a unique R by the confluent hypergeometric functions 
that describe the general solutions of the SCHA problem. 

We note that free IHO energy levels show the well-known "(2z^ + degeneracy with the equidistant eigenvalues 
given by {2iy + £- l)huj, = 1, 2, 3, • • • , for a given £. Such a degeneracy is removed under the confined conditions. 
As E(i',£) > at all R, the critical radius is absent. The incidental degeneracy observed in the case of SCIHO is 
qualitatively similar to that of the SCHA. However, the behavior of the two confined states at a common radius of 
confinement is found to be different [TOl [15]. In particular, for the SCIHO the pairs of the confined states defined 
hy {v = £ + and {v — £ + 2,£ + 2) at the common R = ^ {2£ + 3)/2 a.u., display for all i/, a constant energy 
separation of exactly 2 harmonic-oscillator units, 2huj , with the state of higher £ corresponding to lower energy. The 
choice of R is qualitatively similar to that in the case of SCHA, namely, it is the location of the radial node in the 
{ly ^ £+!,£) state which is the first excited state corresponding to a given £ for the free IHO. It is interesting to note 
that the two confined states at the common R with A£ = 2, considered above, contain different numbers of radial 
nodes. 

With this background, we shall now consider the spherically confined potential defined in Eq.Q. The paper is 
organized as follows. In section 2, the scaling properties and monotonicity of the eigenspectrum generated by the 
potential V(r), as a function of the parameters of the potential, are derived. Analytic energy bounds, derived by 
the envelope method, are reported in section 3: these are found to be useful in guiding the search for very accurate 
values by numerical methods. In sections 4 and 5, we use the asymptotic iteration method (AIM) to study how the 
eigenvalues depend on the potential parameters {a,b,R}, repectively for the free system (R — 00), and for finite R. 
In each of these sections, the results obtained are of two types: exact analytic results that are valid when certain 
parametric constraints are satisfied, and accurate numerical values for arbitrary sets of potential parameters. In 
section 6 we adjoin some more numerical data, obtained by the generalized pseudo-spectral (GPS) Legendre method, 
and present a detailed analysis of the spectral characteristics of the system and their experimental significance. 



II. SCALING AND MONOTONICITY 



Since the potential and the confining box are spherically symmetric, we may write the energy eigenf unctions in the 
form 

vE'(r) = ^rr(0,0), ^(0)=0, (2) 

where r € 3?^ and r = |r|. For finite box sizes R we also require ipiR) = 0. In terms of the atomic units used, 
each discrete eigenvalue depends on three parameters. We shall express this by writing E^£ = E(a,b, R). If we now 
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introduce a scale factor (dilation) cr > into the terms of the Hamiltonian, so that r — ar, then, after multiplying 
the eigenequation Hip = Etp through by cr^, we may derive the general scaling law 

E{a, b, R) = a-^E{aa, a% R/a), a > 0. (3) 

For example, the particular choices a — a~^, a — b^i, and a = R, then yield, respectively, the special scaling laws 

E{a,b,R) = a^E{l,ba-^,aR) = b^ E{ab-i ,l,bi R) = R-^E{aR,bR*,l). (4) 

Thus it would be sufficient to consider just two spectral parameters. 

The eigenvalues E^^i — E{a,b,R) are monotonic in each parameter. For a and b, this is a direct consequence of the 
monotonicity of the potential V in these parameters. Indeed, since dV/da = — 1/r < and dV/db = > 0, it follows 
that 

^^%^<0 and ^^(f'^)>0. (5) 
da ob 

The monotonicity with respect to the box size R may be proved by a variational argument. We shall show in 



section ( III I that the Hamiltonian H is bounded below. The eigenvalues of H may therefore be characterized varia- 
tionally. Let us consider two box sizes, i?i < R2 and an angular momentum subspace labelled by a fixed We extend 
the domains of the wave functions in the finite-dimensional subspace spanned by the first N radial eigenfunctions 
for R = Ri so that the new space W may be used to study the case i? = i?2- We do this by defining the extended 
eigenfunctions so that ipi (r) = for i?i < r < i?2 . We now look at H inW with box size R2 . The minima of the energy 
matrix [{ipi, Hipj)] are the exact eigenvalues for Ri and, by the Rayleigh-Ritz principle, these values are one-by-one 
upper bounds to the eigenvalues for i?2. Thus, by formal argument we deduce what is perhaps intuitively clear, that 
the eigenvalues increase as R is decreased, that is to say 

oR 

From a classical point of view, this Heisenberg-uncertainty effect is perhaps counter intuitive: if we try to squeeze the 
electron into the Coulomb well by reducing R, the reverse happens; eventually, the eigenvalues become positive and 
arbitrarily large, and less and less affected by the presence of the Coulomb singularity. 

For some of our results we shall consider the system unconstrained by a spherical box, that is to say R — 00. For these 
cases, we shall write E,ji — E(a,b). If a very special box is now considered, whose size R coincides with any radial 
node of the R = 00 problem, then the two problems share an eigenvalue exactly. This is an example of a very general 
relation which exists between constrained and unconstrained eigensystems, and, indeed, also between two constrained 
systems with different box sizes. 

III. SOME ANALYTICAL ENERGY BOUNDS 

The generalized Heisenberg uncertainty relation may be expressed [16l[T7] as the operator inequality —A > l/(4r^). 
This allows us to construct the following lower energy bound 



E > £ = min 

0<r<R 



1 a 

g^2 J. 



br^ 



(7) 



Provided b > 0, this lower bound is finite for all a. It also obeys the same scaling and monotonicity laws as E itself 

dr ^ dr 



But the bound is weak. For potentials such as V{r) that satisfy -^{r'^^) > 0, Common has shown ^18j for the ground 



state that (—A) > (l/(2r )), but the resulting energy lower bound is still weak. 



For the unconstrained case R — 00, however, envelope methods p!9H23l 125] allow one to construct analytical upper 
and lower energy bounds with general forms similar to ([t]). In this case we shall write E^£ = E[a, b). Upper and lower 
bounds on the eigenvalues are based on the geometrical fact that V{r) is at once a concave function V{r) — g^^\r^) 
of and a convex function V{r) — g^^'>{~l/r) of —1/r. Thus tangents to the g functions are either shifted scaled 
oscillators above V{r), or shifted scaled atoms below V{r). The resulting energy-bound formulas are given by 



mm 

r>0 



±-^+biP,r) 



< E^Aa.b) < min 

r>0 



(8) 
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where (Ref. [M] Eq.(4.4)) 



Pi ^v, P2 = 2iy-{£ 



(9) 



We use the convention of atomic physics in which, even for non-Coulombic central potentials, a principal quantum 
number v is used and defined by v = n + £ +1, where n is the number of nodes in the radial wave function. It is clear 
that the lower energy bound has the Coulombic degeneracies, and the upper bound those of the harmonic oscillator. 
These bounds are very helpful as a guide when we seek very accurate numerical estimates for these eigenvalues. 

Another related estimate is given by the 'sum approximation' [23" which is more accurate than ([8| and is known to 
be a lower energy bound for the bottom Eg^n of each angular-momentum sub space: in terms of the P's we have for 
these states, P2 = v + }^ = Pi + The estimate is given by 



(a, 6) — min 

r>0 



1 

2r^ 



a 
P^ 



biP^rf 



(10) 



This energy formula has the attractive spectral interpolation property that it is exact whenever a or 6 is zero. The 
energy bounds (Is]) and (10 1 obey the same scaling and monotonicity laws is those of E^i{a,b). Because of their 



simplicity they allow one to extract analytical properties of the eigenvalues. For example, we ca n es timate the critical 
oscillator coupling b that will lead to vanishing energy E — We may estimate b by using (|8|) or ( 10 ). We differentiate 



with respect to r, and use the vanishing of this derivative and of E to obtain the following explicit formula for b 

27 
32 



p4 p2 ' 



(11) 



in which Pa and Pb are to be chosen. If Pa = Pi and Pb = P2, then from ( 10 1 we obtain a good general approximation 
for b. We can also obtain bounds on b. Since E{a, b) is a monotone increasing function of 6, we can state the nature 
of the bounds on b given by formula (111: {i) ii Pa — Pb = Pi = i', the formula yields an upper bound; (ii) if 
Pa — Pb — P2 — 2i/ — (€ + 5 ) , then it is a lower bound; (iii) if z/ = ^ + 1 and Pa — ly and Pb = + ^ , then the formula 
yields a lower bound. We shall state this last result explicitly: for the bottom of each angular-momentum subspace, 
where v — £ + 1, the critical oscillator coupling b yielding i? = is bounded by 



b > 



27 

32 j i.4(^+ 1)2- 



(12) 



IV. EXACT SOLUTIONS FOR THE POTENTIAL V{r) 



The radial three-dimensional Schrodinger equation for the Coulomb plus harmonic-oscillator potential, expressed in 
atomic units, is given by 



1 d'^iPir 

2 dr^ 



l{l + l) 
2r2 



ip{r) ^ E^{r), 0<r<oo, 6>0, aS 



(13) 



where 1(1 + 1) represents the eigenvalue of the square of the angular-momentum operator L^. Note that for a — 0, the 
potential V{r) — —a/r + br^ corresponds to the pure harmonic oscillator potential, while for a > 0, it is a sum of two 
potentials, the attractive Coulomb term —a/r plus the harmonic-oscillator potential br^. For a < 0, the potential V{r) 
corresponds to the sum of two potentials, the repulsive Coulomb potential \a\/r plus a harmonic-oscillator potential 
br\ 

Since the harmonic oscillator potential dominates at large r, this suggests the following Ansatz for the wave function: 

V>(r) = r'+i exp{-ar^)f{r), (14) 



where a is a positive parameter to be determined. Substituting this wave function into Schrodinger's equation ( |13[ ), 
we obtain the following second-order differential equation for f{r): 



rf"{r) + {-Aar^ +21 + 2)/'(r) + ((-25 + 4:a^)r^ + {-ial + 2E- 6a)r + 2a)/(r) = 0, 



(15) 
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which suggest the value a = \/b/2. With this value of a, Eq.(15) is reduced to 
r/"(r) + (-2r2\/26 + 2(; + 1)) /'(r) + [(2^ - (21 + 3)V2b 



2a 



fir) = 0. 



In order to find the polynomials solutions /(r) = X]fc=o ^fe^*^ of this equation, we rely on the following theorem 
Theorem 5) that characterizes the polynomial solutions of a class of differential equations given by 

(03,02;^ + as.ix^ + 03,22; + 03,3) y" + {a2fiX^ + 02,12; + 02,2) y' ~ {nfiX + ri,i) y = 0, 
where 03,,;, i = 0,1, 2, 3, a2,j,j = 0, 1, 2 and ti ^, fc = 0, 1 are arbitrary constant parameters. 
Theorem 1. The second-order linear differential equation | [_?7[ ) has a polynomial solution of degree n if 

Tifl = n{n - 1)03,0 + na2fi, n = 0, 1, 2, . . . , 
along with the vanishing of (n + 1) x (n + \)- determinant A„4_i given by 

l3o ai rii 
71 Pi a-2 772 

72 /?2 ^3 '73 

7n-2 fin-2 "n-l Vn-l 

7n-l an 

In I3n 

where its entries are expressed in terms of the parameters of Eq. by 

Pn = - - 1)03,1 + 02,1) 

a„ = -n({n - 1)03,2 + 02,2) 
In = n,o - (n - l)((n - 2)a3,o + 02,0) 

Vn = -"-("- + 1)03,3 



(16) 



(17) 



(18) 



(19) 



Here, ti.q is fixed by Eq.(18) for a given value of n; the degree of the polynomial solution. 



Consequently, for the polynomial solutions of Eq.(16), we must have, by means of Eq.(18), that 



Eni = in + l+-)V2b 



(20) 



and the conditions on the potential parameters are determined by the vanishing of the tri-diagonal determinant with 
entries 

= -2a 
an = —n{n + 2/ + 1) 
7„ = 2{n-k- l)V2b 
Vn = 

namely, the vanishing of the (n + 1) x (n + l)-tridiagonal determinant 



(21) 



Aji+i — 



-2a -(2 + 2/) 
-2kV2b -2a 

2(1 - k)^/2h 



-2(3 + 20 
-2a 



-3(4 + 20 



2(n-3-fc)V26 -2a -(n-l)(n + 20 

2{n-2-k)V2b -2a -n{n + 2l + l) 

2{n-k-l)V2b -2a 
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For n = we have, for the purely harmonic oscihator a = 0, the exact energy 



Em = (/ + -)V26 (22) 

which gives the ground-state /o(a;) = 1 in each subspace labeUed by the angular momentum quantum number I. 
For n = 1, the determinant A2 — forces the potential parameters a and b to satisfy the equality 



2b{l + 1) = 



(23) 



with a necessary condition for the eigenenergy 



Ell = {l + -W2b. 



(24) 



The condition (23) gives two possibilities for the wavefunction solution. First, for a = — y v26(Z + 1), i.e. with 
repulsive Coulomb term, we have a ground-state (no-node) eigenfunction given by 



(25) 



while for a = a/ V2b{l + 1), i.e. an attractive Coulomb term, we have a first-excited state (one-node): 



tpi{r) — r'"'^^ exp I — 



f - 



26 



(26) 



In table (|T|, we report the first few exact solutions along with the conditions on the potential parameters. Note, the 
subscripts on the polynomial solutions fi{r) refer to the possible number of nodes n in the wave function. 



TABLE I: Conditions for Exact Solutions, here Eni = {n + I + |)v2fo 

Mr) = 1 
a = 



1 /„<o,„.o(r) = l-^, 



a-" - V26(; + 1) = 

a>0,Ti=2 



a(a^ - V2b(5 + 4/)) = 

o f -I g ^ I a^-3^/2b(l + l) 2 l a(a^-y26(7i+9)) 3 

^l<oZ=o l^ ' i+l' {!+l)(2i+3) 3 (I + l){i+2)(2!+3) ' 



a"* - 5\/26(3 + 2l)a^ + 186(2 + + l)=Q 



For arbitrary values of a and b that do not satisfy the conditions (20 1 and (21 1, we may use the asymptotic iteration 



method [27 that can be summarized by the following theorem (for details, see [27], section V Theorem 1, and |28j . 
equations (2.13)-(2.14)): 

Theorem 3: Given Aq = Ao(a;) and sq = sq{x) in C°° , the differential equation 

= Xo{x)y' + so{x)y (27) 
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has a general solution 



I X \ X / t 

y = expi- J a{t)dt j C2 + Ci J exp i J (Ao(r) + 2a(T)) dr | dt 



if for some n > 



where, for n> 1, 



An A„_i 



a{x), or Snix) ^ XnSn-l - Xn-lSn ^ 0, 



Sn — Sn— 1 SoAn- 



(28) 



(29) 



(30) 



Thus, for Eq.(16l, with Ao(r) and SQ{r) given by 

Ao(r) = -i (-2r2x/25 + 2(Z + l)) 



(31) 



so(r) = -i 2£; - (2? + 3)\/26 r + 2a 



the asymptotic iteration sequence A„(x) and s„(x) can be calculated iteratively using (30 1. The energy eigenvalues 
E = Eni of Eq.(16) can be obtained as roots of the termination condition (29). According to the asymptotic iteration 



method (AIM), in particular the study of Brodie et al unless the differential equation is exactly solvable, the 
termination condition ( |29[ ) produces for each iteration an expression that depends on both r and E (for given values of 
the parameters a, h and I). In such a case, one faces the problem of finding the best possible starting value r — r^ that 
stabilizes the AIM process [29]. For our problem, we find that the starting value of rg = 4 is sufficient to utilize AIM 
without much worry about the best possible value of rg. For small values of a, where the wavefunction is spread out, 
we may increase rp > 4. In Table [nj we report our numerical results, using AIM, for energies Eni for the attractive 
(a = 1) and repulsive (a = —I) Coulomb term plus the harmonic-oscillator potential. The numerical computations 
in the present work were done using Maple version 13 running on an IBM architecture personal computer where 
we used a high-precision environment. In order to accelerate our computation we have written our own code for 
root-finding algorithm using a bisection method, instead of using the default procedure 'Solve' of Maple 13. The 
numerical results reported in Table [TT| are accurate to the number of decimals reported. The subscript N refers to the 
number of iterations used by AIM. 



TABLE II; Eigenvalues Eni for V{r) = — a/r -I- br^ , where h = 0.5, a — ±1 and different n and I. The subscript A*' refer to the 
number of iteration used by AIM. 



6 = 0.5 


a = 1 


a = -l 


n 


/ 


Eni 


n 


/ 


Eni 


1 





2.500 000 000 000 000 000N^3,exact 








2.500 000 000 000 000 000N=3,exact 




1 


3.801 929 609 626 278 046iv=80 




1 


3.219 314 119 830 611 360iv=74 




2 


4.930 673 420 047 524 772jv=72 




2 


4.087 227 795 734 562 981iv=67 




3 


6.006 537 298 710 828 780Ar=65 




3 


5.007 681 882 732 318 957iv=6i 




4 


7.058 140 776 824 529 475iv=60 




4 


5.953 327 675 284 371 524jv=56 








0.179 668 484 653 553 873iv=97 


1 





4.380 233 836 413 610 273iv=97 


1 




2.500 000 000 000 000 OOOn^s, Exact 


2 




6.301 066 353 339 463 595iv=67 


2 




4.631 952 408 873 053 214jv=72 


3 




8.243 517 978 923 477 298iv=67 


3 




6.712 595 725 661 429 760iv=7o 


4 




10.199 062 810 923 865 963]v=65 


4 




8.769 519 600 328 899 714jv=69 


5 




12.163 259 523 048 320 928]v=64 
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V. EXACT SOLUTIONS FOR THE SPHERICALLY CONFINED V{r) 



In this section, we consider the confined case of Coulomb and liarmonic oscillator system as described by the radial 
Schrodinger equation (in atomic units) 



2 dr2 



2r2 



--+br' 
r 



i/,(r) = E^Pir), 0<r <R, 6 > 0. 



(32) 



where / = 0, 1, . . . is the angular-momentum quantum number and '0(0) — V'(^) — 0- Here, again, the parameter a is 
allowed to in M. Intuitively, we may assume the following ansatz for the wave function 

iP{r) = r^+\R - r) exp{-ar^ - I3r)f{r), (33) 

where a and /? are parameters to be determine, and R is the radius of confinement. The R — r factor ensures that the 
wave function will become zero at r = i?. Direct substitution of Eq.(33) into Eq.(32) yields the following second-order 
linear differential equation for /(r): 

'/ + 1 1 



nr) ^ -2 



1 



R — r J 

(25 - Aa^y^ + {ARa^ - Ajia - 2Rh)r^ + {ARaf3 - 2E - (5^ + Ala + Way 



r{R - r) 

{Rp^ - 6Ra + 2RE - ARla + 21(3 + A^ - 2a)r - 2{l + 1) + 2Ra - 2i?/3(Z -t- 1 



fir) 



(34) 



Clearly, from this equation, we have a — \/b/2 and /? = 0, which reveals the domination of the harmonic oscillator 
term even in the confined case. Consequently, for /(r), we have 



R-r 



2b r fir) 



r{R - r) 



i-2E + {21 + 5)V2by + {-3RV2b + 2RE - 2RlV2b - 2a)r - 2{l + 1) + 2Ra 



fir) (35) 



Although, equation (35 1 still does not lie within the framework of Theorem 1, we may make use of the following result 
(HHj, Theorem 6) 

Theorem 4. A necessary condition for the second-order linear differential equation 

/k+2 \ /fc+l \ / k \ 



E 

Vi=0 / \j=0 / \i=0 

to have a polynomial solution of degree n is 

Tk,o = ri{n " l)afe+2,o + "■a/c+i.o, A; = 0, 1, 2, 



TkAX 



k—i 



y = 



(36) 



(37) 



Thus for Eq.(35), or, more explicitly, the differential equation 



r{R - r)f"(r) + 2 ((Z + l)(i? - r) - r - V2b r^R + V2b r^)) /'(r) 

-I- {-2E + {21 + 5)V2b)r^ + {~3RV2b + 2RE - 2RlV2b - 2a)r + 2Ra - 2{l + 1] 



to have polynomial solutions of the form /„(?■) = X]fc=o '^k^'^, it is necessary that 



E„i = {n + l+-)V2b. 



fir) = (38) 



(39) 



This is an important formula for Eni that can facilitate greatly our computations based on AIM. We note, first, using 
Eq.(39l that Eq.(35l can be reduced to 



f::ir) 



i+i 



1 



R-r 



y26r /;(r) 



1 

r{R-r) 



- 2nV2b r^ + (2i?V2&(n + 1) - 2a)r + 2Ra - 2{l + 1) 



fair). 



(40) 
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It is then clear from equation (401 that, for n = 0, we have 

r Eoi ^ (l + l)V2b, /o(r) = l, 



[ V'o;W =r'+i(i?-r)exp(-^|r2). 
if the parameters a, b and the radius of confinement R are related by 

aR^l + 1, = + l)%/2&. 



(41) 



(42) 



The wavefunction given by (41 1 represent the ground-state eigenfunction is each subspace labeled by the angular 
momentum quantum number I. 

• For n = 1, we, easily, find that 

< 

only if the parameters a, b and R are related by 



(43) 



^ I _ 1+1 



^ ° ^ 2 I J 



yR R^ I ' 

« = if (2^ + I ± 1^41^5) . 
Or, more explicitly, for a and b expressed in terms of the radius of confinement R, as 

' a± = ;^(2/ + |± iV'¥T5), 



(44) 



E = ^(7 + 2Z)\/26. 



(45) 



From (431, for a > 0, it is clear that either a < [l + or a > (^ + since, for the case of a = (/ + l)/i?, 

we have 6 = 0, which is not acceptable from the structure of our wave function (33) where 6 > 0. We further 



note from (43) that for r < i? to have one node within (0,i?), it is necessary that R> 2(1 + l)/a > (^ + l)/a. 
For example, if a = 1,Z = 0, then for a one-node state within (0, i?), it is required that R > 2. Thus, let 
a = 1,1 = 0,R = 5/2 + V5/2, we have from b 1/50 and E^ = 0.700 000 000 000 000. Note further, 

if a = 1, ; = but i? = 5/2 - ^5/2, although 1/R < (1 + l)/a, still we do not have any node that lies within 
(0, 5/2 — a/5/2), since in this case r = R = 5/2 — \/5/2. Thus we have, in this case, a node-less wave function 
/o(r) with Eoo = 0.699 999 999 995 412 275. This explains the subscript /o,i in (|43). 



Further, for n = 2, we can show that 



' E2i = l{9 + 2l)V2b, /o,i,2(r) = 1 + (ij - i^t) ^ " ( 
^ V2/(r) = r'+i(i? - r) cxp(-^^r') (l + (l - ^) r - 



i3V2bil+l)~a'^)R^ + {aR-l-l){2l+3) 
fl2(; + l)(2i+3) 



ft2(; + l)(2i+3) 



(46) 



only if a, b and R are related by the two-implicit expressions: 



and 



R-^a-^ - ZR\l + 2)a^ - R{RW2b{7l + 9) - 3(Z + 2){2l + 3))a - 3(/ + 2){l + l){2l + 3 - 3i?V26) 



R^a^ - R{R'^V2b + 21 + 3)a2 + (l + 1){21 + 3 - m^V2b)a + 6bR^{l + 1) = 0. 
We now consider a few examples of these results: 



(47) 
(48) 
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If a = IJ = 0, then (47) and (48 1 yields V2b = 0.76025880213480504582, i? = 2.0843217092058454961, 



and we have the following solution: 

Eio = 3.421 164 609 606 622 706 2, 

i/;io(r) = r(2.084 321 709 205 845 496 1 - r)exp(-0.380 129 401 067 402 522 91 r^) 
X (l - 0.520 227 613 816 384 707 07 r - 0.676 516 312 440 766 915 00 r^) . 



(49) 



If V2b = 2,/ = 0, then (47) and (48) yields a 



and we have the following solution: 



-1.6219380762368883824, i? = 1.0232416568868508038 



Eqo = 9.000 000 000 000 000, 



Vjoo(r) = r(1.023 241 656 886 850 803 8 - r) exp(-r2) 

X (1 + 2.599 224 324 572 826 833 6 r + 1.417 080 563 127 630 983 r^) 



(50) 



• For n = 3, we have 

' E31 = i(ll + 20^26 



1 {-R^a^+3R.^{l+2)a^+B^F/.^^/2b{10l + 13)-3(l+2){2l+3))a-3il+2){l + l){-2l+4R^^/2b-3)) 3 
+ 3 m(l+l){l+2){2l+3) 



_ R(4R.V2b{l + l)-a{Ra-3-2l))-{2l+3){l+l) 2 
i?2(2;+3)(i + l) ^ 



(51) 



where a, b and R are related by 

R^a^ - R^{R^V2b + 3(/ + 2))a^ - R{R^V2b{l0l + 13) - 3{l + 2){2l + 3))a 



and 



26(13 + 10Z)i?'* + 12V26(; + 1)(/ + 2)R^ - 3(2Z + 3)(Z + 2){l + 1) = 



R'^a^ - 2R^{21 + 5)a3 - R^{R^V2b{m + 25) - 6(2/ + 5)(/ + 2))a^ + 2R(2l + 5){R^V2b{10l + 13) 
- 3(/ + 2)(2/ + 3))a + 6(/ + 1)(/ + 2)((2/ + 5)(2/ + 3) + 4,V2bR^ {R^ V2h - 5 - 2/)) = 0. 



(52) 
(53) 



(54) 



Similar results can be obtain for higher n (the degree of the polynomial solutions) . It is necessary to note that the 
conditions reported here are for the mixed potential V{r) = —ajr + 6r^, where a 7^ 0, 6 7^ (neither coefficient is 
zero) . 

For the arbitrary values of a, b and i?, not necessarily satisfying the above conditions, we still apply AIM directly 
to compute the eigenvalues. Similarly to the un-confined case, we start with 



Ao(r) 
So(f) ^ 



r{R-r) 



2b r) , 

[-2E + (2/ + 5) V25)r2 + {-3RV2b + 2RE ~ 2RlV2b - 2a)r + 2Ra - 2{l + 1) 



(55) 



The AIM sequence A„(x) and Sn{x) can be calculated iteratively using (30). The energy eigenvalues E = Eni of 



Eq.(38l are obtained as roots of the termination condition (31). Since the differential equation (|38|) has two regular 



singular points at r = and r = R, our initial value of ro can be chosen to be an arbitrary value in (0, R). In table III 



we reported the eigenvalues computed using AIM for a fixed radius of confinement R — 1 with tq = 0.5 as an initial 
value to seed the AIM process. In general, the computation of the eigenvalues are fast as illustrated by the small 
number of iteration N in Tables |III[ |IV| an d [V] The same procedure can be applied to compute the eigenvalues for 
arbitrary values of a, b and R. In Table IV we have fixed a, b and allowed R to vary, then we fixed 6, R and allowed a 



to vary. In Table [Vj we fixed a, R and varied b. Our numerical results in these tables confirm our earlier monotonicity 
formulas reported in section II. 
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TABLE III: Eigenvalues Eni for V{r) = -a/r + br^ , r G (0,i?), where h = 0.5, a = ±1, i? = 1 and diflterent n and I. The 
subscript TV refers to the number of iteration used by AIM. 



a= 1,6 = 0.5, i?= 1 


n 


I 




n 




Enl 








2.500 000 000 000 000 QQ0N=3,Exact 








2.500 000 000 000 000 OOOjv^s.B^act 




1 


8.404 448 391 842 929 575jv^24 


1 




16.733 064 961 893 308 967iv=25 




2 


15.183 570 193 031 143 001jv=23 


2 




41.029 002 263 262 675 364jv=33 




3 


23.137 256 709 545 767 885iv=24 


3 




75.297 038 665 283 580 892jv=4o 




4 


32.295 207 272 878 341 541jv=27 


4 




119.493 804 921 354 632 859iv=47 


a = -1,6 = 0.5, i? = 1 


n 


Z 


E„i 


n 


/ 


Enl 








7.427 602 986 235 605 737jv=26 








7.427 602 986 235 605 737jv=26 




1 


12.118 629 877 542 593 085iv=24 


1 




22.954 866 627 528 634 394jv=27 




2 


18.456 796 172 766 948 526jv=23 


2 




48.054 781 032 847 609 425iv=36 




3 


26.173 002 039 626 403 748iv=25 


3 




82.897 495 765 909 946 966iv=4i 




4 


35.179 533 437 869 611 594jv=28 


4 




127.540 759 830 804 826 131iv=48 



TABLE IV; Eigenvalues Eoa for V{r) = -a/r + 6r^ r G (0, 7?), where we fixed 6 = 0.5 and we allowed a and R to vary. The 
subscript TV refers to the number of iteration used by AIM. 









6 = 


0.5 








a 


R 


Eoo 


R 


a 


-Boo 


1 


0.1 


468.994 438 340 395 273 843iv 


=26 


1 




10 


24.446 394 090 129 924 468iv= 


=25 




0.5 


14.781 525 455 450 240 772jv. 


= 19 






-5 


15.581 919 590 917 726 881]v= 


=25 




1 


2.500 000 000 000 000 000iv= 


Exact 






-1 


7.427 602 986 235 605 737jv= 


26 




2 


0.281 457 639 408 567 801iv= 


44 









5.075 582 015 226 783 066jv= 


26 




3 


0.180 768 103 642 728 017iv= 


66 






1 


2.500 000 000 000 000 000jv= 


Z, Exact 




4 


0.179 669 842 444 710 526iv= 


80 






5 


- 12.356 931 301 584 560 963jv 


=35 




5 


0.179 668 484 856 687 713iv= 


82 






10 


- 49.984 937 021 677 890 425jv 


=43 



TABLE V: Eigenvalues £00 for V{r) = -~a/r + br^, r G (0, R), where we fixed a = R = 1 and allowed b to vary. The subscript 
TV refer to the number of iterations used by AIM. 



7? = l,a = 1 


b 


Eoo 


0.1 


2.399 281 395 696 719 214iv= 


=22 


0.2 


2.424 527 479 482 894 839]v= 


=22 


0.5 


2.500 000 000 000 000 000]v= 


=3 


1.0 


2.624 907 458 899 526 414]v= 


=31 


2.0 


2.871 465 192 314 860 746]v= 


=35 


5.0 


3.585 958 081 033 459 432]v= 


=41 


10.0 


4.698 782 960 476 752 179jv= 


=47 



VI. SPECTRAL CHARACTERISTICS 

In this section we shall discuss the spectral characteristics associated with the crossings of the energy levels. We have 
employed the generalized pseudo-spectral (GPS) Legendre method with mapping, which is a fast algorithm that has 
been tested extensively and shown to yield the eigenvalues with an accuracy of twelve digits after the decimal. A 
more detailed account, with several applications of GPS, can be found in [SUMS] and the references therein. 

In the present work, we have also verified the accuracy of these results, in a few selected cases, by using AIM. We 
shall first consider the case defined by i? — t- 00, a — 1 and variable b, under which the potential given by Eq.(l) can 
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be regarded as representing the hydrogen atom confined by a soft harmonic oscillator potential. Starting from the 
free HA (6 = 0), the effect of finite b is to remove the accidental degeneracy and raise the energy levels such that 
Eliy, t) > E{i^,£+1) (see Ref. [SS])- As the starting Elv, £) < given by the HA spectrum, there exists a critical value 
of b — be, corresponding to each level, defined by the condition E(y,t) — 0. The numerical values of be are found 
to be rather small, except for the ground state, indicating that a weak confinement due to the harmonic potential is 
sufficient to realize the condition that E{v,€) > for all b > be- In the usual spectroscopic notations the levels 

(Is2s2p3s3p3d4s4p4(i4/) 

are defined by the be values given respectively by 

be = (0.32533, 0.004831, 0.00771, 0.00042, 0.00051, 0.00079, 0.00007, 0.00008, 0.00010, 0.00015). 

In Fig. 1, we have displayed the passing of the energy levels corresponding to AsApAdAf states through E — a.t 
be- We know tht the eigenspectrum of free HA is indeed very sensitive to the harmonic confinement since it is 
found numerically that at 6 = 0.000001, the eigenvalues are already positive, corresponding to the states given by 
7p, 7d, 8d, .., 7/, 8/, .., 7(;... The crossings of energy levels can be gauged by the change of ordering from the hydrogen- 
like 

{ls2p2sMSp3s4:f4:d5g'ip4:s5f5d6g5p5s6f6d7g6p6s7f7d8g7p8f8d9g9f ...) 

to 

^ {Is2p2s3d3p4:f3s4d5g4:p5f4:s5d6g5p6f5s6d7g6p7f6s7d8g7p8f8d9g9f . . . ) 

as the parameters of the potential change along (a = 1, 6 = 0) — > (a = 1, 5 = 0.001) — > (a = 1,6 = 0.5). It follows 
that the (3s, 4/), (4s, 5/) . . . levels defined by (i^, £) and + 1,£ + 3) cross at a certain b. 




-0.04 J 
0.00001 



FIG. 1: The critical b, denoted as be in the text at which E{i/, £) — are shown for the 4s, 4p, id, 4/ states. The large value of 
R — 100 corresponds essentially to the free state of the potential in Eq.(l) with a — 1. 

In Fig. 2, we have displayed this behavior corresponding to a = 1. This spectral characteristic is similar to that 
found earlier [37] for the case of the soft Coulomb potential. Further, the eigenvalue [v = ^,£ — A) is found to cross 
{v -l,£- 4),{v -\,£-i){v,£- 4),{v,£~'2,){v,£-2),{v,£-l) asb changes from to 0.5. 
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b 



FIG. 2: The crossings of levels as a function of fe, corresponding to the free state of the potential in Eq.(l) with a — 1. The 
levels defined by (Vji) and (i/ + 1,£ + 3) are shown. 

Next, we consider the new spectral characteristics introduced when, in addition to the harmonic-oscillator potential 
term, a second confining feature consisting of an impenetrable sphere of finite radius R is introduced. Such a potential 
factor further raises the energy levels as R is diminished, — > 0. As a consequence the be values get smaller. This is 
depicted in Fig. 3 for the 4s and Ap states at two different values of R of 100 and 30 a.u., respectively. The former 
corresponds to the case i? -> oo, i.e. just the potential in Eq.(l). Varying R under fixed a yields a different level 
ordering, depending upon the value of b, as this situation corresponds to two specifically chosen confinement features 
imposed on the hydrogen-like potential at each point. To illustrate this, we consider the case defined by a = 1, & = 0.5 
and variable R. Our calculations suggest that the ordering of levels changes from 

{ls2p2s3d3p4fSs4:d5gAp5fAs5d6g5p6f5s6d7g6p7f6s7d8g7pSf8d9g9fl0g ...) 

to 

^ (Is2j93d2s4/3p5g4d3s5/4p6g5d4s6/5p7g6d5s7/6p8.g7d6s8/7p958d9/10g . . . ) 

as R changes from oo 0. The crossings of levels are now observed between the state {v, £) and {v — 1,£ + 2). 
In Fig. 4, we have shown this feature corresponding to the confined {3s, M) and (3s, 4/) states. Additionally, the 5g 
level is found to fall below 4d and 3s levels, successively, as R decreases. It is evident that the imposition of a double 
confinement effect, mediated through the combination of and the boundary at R leads to the crossings among a 
wider set of the states of the hydrogen-like atom, not observed in the separate singly confined situations. A possible 
experimental system of embedded atom inside zeolite, fuUerine, or liquid helium droplets under very strong laser fields 
could be modelled using the doubly confined Coulomb potential as described in this work. 
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E=0 at Critical b witti a=1, R=30 and R=100 




-4s(100) 
-4s(30) 
-4p(100) 
-4p(3Q) 



FIG. 3: The critical b, denoted as be in the text at which E{iy, £) = are shown for the 4s, 4p states. The value of be decreases 
as R decreases: specifically, the essentially free state of the potential in Eq.(l) with a = 1 at i? = 100 is confined to a smaller 
value of i? = 30. The numbers inside brackets denote R. 



Crossing of energy levels as a function of R 
[a=1,b=0.5] 




FIG. 4: The crossings of levels as R is changed as the potential in Eq.(l) is defined by the values a = 1 b — 0.5. Crossings are 
observed between the state {ly, £) and {i/ — 1,£ + 2) as shown by the (3s, 4ci) and (3s, 4/) levels. The 5g level is shown to cross 
4d and 3s as R decreases. 
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VII. CONCLUSION 



In this study we first consider a very elementary model for an atom, namely a single particle which moves in a central 
Coulomb potential — a/r and obeys quantum mechanics. We then adjoin two confining features: soft confinement by 
means of an attractive oscillator term br^, and hard confinement produced by containment inside an impenetrable 
spherical cavity of radius R. The paper reports on the effects of the confinement parameters {b, R} on the original 
Coulomb spectrum which, of course, is given in atomic units by i;^ = —a^ /{2u'^). By a combination of analytical 
and numerical techniques, we are able to make considerable progress in analyzing the spectral characteristics of this 
confined atomic model. In future work we plan to undertake a similar study in which the pure Coulomb term is 
replaced by a more physically interesting screened-Coulomb potential, or a soft-core potential such as —a/{r + j3). 
The purpose of this work is to look at model problems that contain physically interesting features but are still simple 
enough to yield to analytical as well as purely numerical analysis. 
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